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ABSTRACT 

Context. Baryon Acoustic Oscillations (BAOs) are oscillatory features in the galaxy power spectrum which are used 
as a standard rod to measure the cosmological expansion. These have been studied in Cartesian space (Fourier or real 
space) or in spherical harmonic space in thin shells. 

Aims. Future wide-field surveys will cover both wide and deep regions of the sky and thus require a simultaneous 
treatment of the spherical sky and of an extended radial coverage. The Spherical Fourier-Bessel (SFB) decomposition 
is a natural basis for the analysis of fields in this geometry and facilitates the combination of BAO surveys with other 
cosmological probes readily described in this basis. In this paper, we present a new way to analyse BAOs by studying 
the BAO wiggles from the SFB power spectrum. 

Methods. In SFB space, the power spectrum generally has both a radial (fe) and tangential {£) dependence and so do 
the BAOs. In the deep survey limit and ignoring evolution, the SFB power spectrum is purely radial and reduces to 
the Cartesian Fourier power spectrum. In the opposite limit of a thin shell, all the information is contained in the 
tangential modes described by the 2D spherical harmonic power spectrum. 

Results. We find that the radialisation of the SFB power spectrum is still a good approximation even when considering 
an evolving and biased galaxy field with a finite selection function. This effect can be observed by all-sky surveys with 
depths comparable to current surveys. We also find that the BAOs radialise more rapidly than the full SFB power 
spectrum. 

Conclusions. Our results suggest the first peak of the BAOs in SFB space becomes radial out to f ~ 10 for all-sky 
surveys with the same depth as SDSS or 2dF, and out to ^ ~ 70 for an all-sky stage IV survey. Subsequent BAO peaks 
will also become radial, but for shallow surveys these may be in the non-linear regime. For modes that have become 
radial, measurements at different ^'s are useful in practice to reduce measurement errors. 

Key words. Cosmology, Baryon Acoustic Oscillations, Wide-field surveys, Statistics, Spherical Fourier-Bessel. 



1. Introduction 

The study of Large Scale Structure (LSS) with galaxy 
surveys is a promising tool to study the dark universe 
(Peacock et al. 2006; Albrccht et al. 2006). Baryon Acoustic 
Oscillations (BAOs) are a special feature in the galaxy 
power spectrum present on scales 100/i~^Mpc, which are 
due to oscillations in the coupled baryon-photon fluid be- 
fore recombination (Sunyaev & Zeldovich 1970; Peebles & 
Yu 1970; Eisenstein et al. 2005; Seo & Eisenstcin 2003, 
2007). The BAOs are considered a powerful cosmological 
tool as the BAO scale acts as a standard ruler with which 
to probe cosmic expansion both in the radial and tangential 
directions. 

BAOs were first detected by Eisenstein et al. (2005) in 
SDSS data (Adelman-McCarthy et al. 2008) and later with 
2dF galaxies (CoUcss ct al. 2003; Cole et al. 2005) and fi- 
nally with both surveys (Percival et al. 2007 b) , though oth- 
ers suggest current data cannot currently probe the BAO 
scale sufficiently (Cabre & Gaztafiaga 2011). 

Until now, BAOs have been studied in Cartesian space, 
either in Fourier space (Seo & Eisenstein 2003, 2007), or in 
real space (Eisenstein et al. 2005; Xu ct al. 2010; Slosar ct al. 



2009) , and in 2D spherical harmonic space on thin spherical 
shells (Dolney et al. 2006). These descriptions use different 
information and therefore have different constraining power 
for cosmological parameters (Rassat et al. 2008). 



Future wide-field BAO surveys will, however, cover both 
large and deep areas of the sky, and thus require a si- 
multaneous treatment of the spherical sky geometry and 
of extended radial coverage. The Spherical Fourier-Bessel 
(SFB) decomposition is a natural basis for the analysis of 
fields in this geometry. The SFB analysis is powerful as it 
uses a coordinate system in which the radial selection func- 
tion and physical effects are naturally described. Moreover, 
this description facilitates the combination of BAO sur- 
veys with other cosmological probes which are readily de- 
scribed in the SFB decomposition such as the smooth power 
spectrum and redshift space distortions (Heavens & Taylor 
1995; Fisher et al. 1995b; Percival et al. 2004; Erdogdu (a) 
et al. 2006; Erdogdu (b) et al. 2006), weak lensing (Heavens 
2003; Castro ct al. 2005; Kitching ct al. 2008), and the 
Integrated Sachs-Wolfe effect (ISW, Shapiro et al. 2011). 
Studying BAOs from wide-field surveys with an SFB ex- 
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pansion, is therefore natural both for the geometry consid- 
ered, and for unifying the treatment with the other probes. 

In §2 we first review the general decomposition of a 3D 
field in SFB space and introduce the concept of radiali- 
sation when the field is statistically isotropic and homo- 
geneous. In §3.1, we introduce a new way to consider the 
BAOs, by studying the baryon wiggles in SFB space, in a 
similar way as in Fourier space, while in §3.2 we discuss the 
radialisation of BAOs in SFB space in the context of exist- 
ing and future surveys. In §4 we present our conclusions. 

2. 3D Spherical Fourier-Bessel (SFB) Expansion 

2.1. Expansion of an homogeneous and isotropic field 

Let us consider a field /(r) at time t, where r = {r,6,Lp) 
in spherical polar coordinates. In practice, this field may 
represent the galaxy or mass density (or overdensity) in the 
universe. In a flat geometry, the field can be decomposed in 
the 3D SFB basis set, which is complete and orthonormal, 
as 



fir) 



dky^Jim{k)kje{kr)Ye„^{6,(p), 



(1) 



im 



where je{x) are spherical Bessel functions of the first kind, 
Y£m{0,ip) are spherical harmonics, i and m are multipole 
moments and k is the wavenumber. The inverse relation is 



fhn{k) 



A^v f[v)kn{kr)Yl,^{e,^), 



(2) 



where we use the same conventions as Leistedt et al. (2011) 
and Castro et al. (2005) (see also Heavens & Taylor 1995; 
Fisher et al. 1995b, who use a different convention and basis 
set). 

This decomposition can be viewed as the spherical po- 
lar coordinate analogue to the Fourier decomposition in 
Cartesian coordinates given by 

/W = /d'fc /»e'''•^ (3) 



(4) 



The 3D SFB power spectrum Ci{k) of the field fir) is 
given by the the 2-point function of the SFB coefficients 
fe.m{k) which can be written as 

{ftm{k)ftra'{k')) = C,{k)5{k ~ k')Su'S^ra', (5) 

when the field is Statistically Isotropic and Homogeneous 
(hereafter, the SIH condition). Similarly and in the same 
condition, the Fourier power spectrum P{k) is implicitly 
defined by 



(27r)3 J 
/(k) = / d'x /(x)e-*'^-. 



imrik')) = {2nfPik)S'ik~k'). 



(6) 



These two power spectra are related by (see for e.g.. 
Fisher et al. 1995b; Castro et al. 2005): 

C,{k)^P{k). (7) 

Thus, the SFB power spectrum Ce(k) is independent of the 
multipole £, and thus only has a radial (k) dependence. This 
remarkable yet recondite fact is only true if the field fulfills 
the SIH condition. 

In the following subsections, we discuss the impact of 
the (partial) violation of this condition in practice and the 
implication for the measurements of BAO in SFB analysis. 



2.2. Finite Surveys 

In practice, a cosmological field, such as the galaxy density 
field, will only be partially observed due the finite survey 
volume. In this case the observed field /°'^^(r) can be de- 
scribed by 



r%v) = 0(r)/(r), 



(8) 



where 0(r) is the radial selection function of the survey, 
and /(r) is assumed to fulfill the SIH condition for now 
(we will address the effect of bias and evolution in §2.3). 
The observed field /°'^''(r) is thus no longer homogeneous 
because of the radial the selection function. There may also 
be a tangential selection function to account for regions of 
missing data, but we assume here that the data is available 
on the full sky. For convenience, the selection function is 
normalised as 



d^r 



V, 



(9) 



where y is a characteristic volume of the survey chosen 
such that (/) — >■ 1 as y — )■ 00, such that 



hm r^^r) - fir). 



(10) 



In this case, the homogeneity condition (SIH) is no 
longer valid in the radial direction and the observed 2-point 
function can be written as 



(11) 



where the observed SFB power spectrum Cg^^{k,k') de- 
pends this time on £, k and k', as opposed to only on £ and 
k as in Equation 5. It can be shown that it can be expressed 

as 



Cr'(fc,/c') 



dfc" k"^P{k")Wi(k, k")Wi{k', k"), (12) 



where P{k) is defined in Equation 6 and the window func- 
tion Wi{k,k') is defined as 

We{k,k')= J dr r^(j){r)kjfXkr)je{k'r). (13) 

In practice, Cf^^{k^k') tends to fall off rapidly away from 
the diagonal k = k' and we will often only compute C^^^{k) 
defined by: 



Cr'(fc) = Cf%k,k) 



(14) 



We next evaluate the window function and the observed 
power spectrum for three special cases for the selection 
function. 



2.2.1. Gaussian selection function 

As a first example, let us consider a Gaussian selection func- 
tion defined as 



<^(r) 



(15) 



where tq is a radius parameter and the normalisation obeys 
Equation 9 (with V = tt^Tq). In this case, the window 
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function can be integrated analytically, and Equation 13 
becomes: 



We{k,k') 



k 



U2 



rlkk' 



(16) 



where Iu{x) is the modified Bessel function of the first kind. 
This analytical form considerably facilitates the evaluation 
of Equation 12 (see §A for a way to numerically calculate 
the above Equation for large arguments of Ii,(x)). 

2.2.2. Radial Limit 

In the case where the radial selection function corresponds 
to full radial coverage (i.e. 4>{r) = l,Vr), the window func- 
tion becomes 



W,{k,k')^^5{k~k') 



(17) 



Inserting this expression into Equation 12, we recover 
Equation 7, namely 



Cf%k, k') = Ct{k)6{k - k') = P{k)5{k - fc'). 



(18) 



meaning that the 3D spherical spectrum is only dependent 
on radial coordinate k, as discussed above. 

As we show in the appendix §A, this limit can also be 
obtained by taking the limit rp — > oo for the Gaussian 
weight function of Equation 15 and is achieved in practice 
when the condition (see Equation A. 4) 



(19) 



is satisfied. We refer to this as the radialisation of the field 
and discuss this more in §3. 

2.2.3. Tangential Limit 

In the other extreme case, the radial selection function cov- 
ers only a thin shell of the field at a distance r» 



") = r^,S{r — r*), 



(20) 



where the normalisation of Equation 9 was chosen with 
V — Airr^. Equation 12 then becomes 

Cf'^ik^k') = kk'j,[kT,)j,{k'r,) 

X / <:\k"k"'P{k")[rlj,{k"u)]\ (21) 



This can be related to the statistics of the 2D projected 
field 

/2D(0,^)= j drr'cj,{r)f{r,e,^), (22) 
whose spherical harmonic decomposition is given by 

CO i 

f^{9,^) = E E fi^YUO,^), (23) 

]Z = [ dnf^ie,ip)Y;^{0,ip), (24) 



and whose 2D spherical harmonic power spectrum is 
defined by 



if em ft m') = Q ^U'S, 



(25) 



This power spectrum is related to the 3D Fourier space 
power spectrum P{k) by: 

= ^ y dkk^Pik) [W^''ik)]\ (26) 

where the 2D window function Wi^{k) is given by 

W^ik) = I drr^(t){r)je{kr). (27) 



Note that due to our choice of normalisation for (/)(r) in 
Equation 9, W^^{k) has units of volume. 

In the case of a thin-shell selection function (Equation 
20), this reduces to 



^™ = IJ dkk'Pik) [rljiikr,)]' 



(28) 



Thus, the SFB power spectrum (Equation 21) is simply 
related to the 2D spherical harmonic power spectrum^ by 



Cf\k,k') = -kk'j,{ku)j,{k'r,)Cf. 



(29) 



The term kk' i^{kr^,)ji{k'r^,) is a geometric factor which 
does not depend on the field /. Thus, in this case, all the 
information in the 3D SFB spectrum Ci{k^ k') is contained 
in the 2D power spectrum C]^ , i.e. depends solely on t. 

2.3. Evolution and bias 

For a realistic galaxy field, the SIH condition will not be 
met due to the radial evolution of the field from the growth 
of structure and redshift dependent bias. 

For a biased and evolving field in the linear regime, the 
galaxy power spectrum P°™'(fc, r) will have an explicit dis- 
tance dependence: 



P''^°\k,r)^h^{r)D^{r)P{k), 



(30) 



where the term h(r) is the linear bias, assumed to be 
scale-independent, D{r) is the growth factor and P{k) = 
P{k, r = 0) is the matter power spectrum at r = z = 0. The 
linear approximation will hold up to a redshift-dependcnt 
scale k„ia,x{z). In the standard cosmological model, we find 
kn,a,^(z = 0) ~ 0.12/iMpc"^ (i.e., before the second BAG 
peak) and k^axiz = 2) ~ 0.25/iMpc~"^. 

For the realistic galaxy field, an approximation for the 
spherical 3D power spectrum for the galaxy field is still 
given by Equation 12, but the window function will incor- 
porate the linear growth and bias functions as 

W^^°\k',k)= [ dr r^(j){r)kD{r)b{r)je{kr)je{k'r). (31) 



^ Equation 29 is conceptually different from Equation B3 in 
Kitching et al. (2011): the former is an exact solution for a thin 
shell, while the latter uses the Limber approximation. 
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Wide and Shallow Survey Wide and Deep Survey 




1 10 1 10 



Fig. 1. Ratio Rf (k) of SFB spectrum with and without the physical effects of baryons in {£, k) space for a wide and 
shallow survey (left, — 100/i~^Mpc) and a wide and deep survey (right, rg = 1400/i~^Mpc) using a Gaussian selection 
function in r. The baryonic wiggles are seen both in the radial (fc) and the tangential (£) directions, but as Tq is increased, 
only radial wiggles - with increased amplitude - persist. This radialisation of the information is due to mode-cancelling 
along the line of sight for deep surveys. 



3. Application to Baryon Acoustic Oscillations 

3.1. SFB Analysis of BAOs 

The matter power spectrum in Equation 12 includes the 
physical effects of baryons, which create oscillations in 
Fourier space (Sunyaev & Zeldovich 1970; Peebles & Yu 
1970; Seo & Eisenstein 2003; Eisenstein ct al. 2005; Seo 
& Eisenstein 2007). These Baryon Acoustic Oscillations 
(BAOs) can be isolated by considering the ratio R^{k) in 
Fourier space: 



R^ik) = 



P'^(fc) 



(32) 



where P^{k) is the galaxy (or matter) power spectrum in- 
cluding the physical effect of BAOs, and P'^°^{k) is broad 
band or 'smooth' part of the galaxy (or matter) power spec- 
trum. In linear theory the growth and bias terms cancel out 
so that Equation 32 is independent of z. These oscillations 
can also be probed in 2D spherical harmonic space (Dolney 
et al. 2006). We are interested here to see if they can simi- 
larly be probed in SFB space as well. 

In SFB space, we consider the ratio Rf{k) given by: 



^ C° • (k) 

^£ (^) /=,obs,nob/7 \ ' 

[k) 



(33) 



where similarly to Fourier space, C'^^'^'^{k) is the diagonal 
SFB power spectrum (Equation 14) including the physical 
effects of baryons and (fc) is the 'smooth' part of 

the SFB power spectrum. These are calculated by using 
P^{k) or P"°^(fc) instead of P{k) in Equation 12 and then 
considering Equation 14. The SFB decomposition suggests 
that the BAOs can manifest themselves in both k and £ 
space. 

In Figure 1, we plot Rf{k) for two wide- field galaxy sur- 
veys with a shallow (left) and deep (right) galaxy Gaussian 
selection functions with tq — 100/i^^Mpc, and 1400ft.^^Mpc 
respectively (see §2.2.1). We chose a fiducial concordance 
cosmology with ilm = 0.25, JIde = 0.75,51b = 0.045, 
Wo = -0.95, Wa = 0.0, h = 0.70, as = 0.80, = 1 and 
no running spectral index. We consider both to be 'realis- 
tic' surveys, i.e. with evolution due to growth taken into 
account and a linear galaxy bias (taken as b{z) = \/l + z), 
using Equation 31 for the window function. 

In the LHS of Figure 1 (narrow survey), the BAOs de- 
pend on both k and £ modes, illustrating that the wiggles 
can be measured simultaneously in the £ and k directions, 
which is the first result of this paper. 

However, in the RHS of Figure 1 (wide survey), the 
BAOs appear to have only a radial (fc) dependence, except 
at very large physical scales {k < 0.01/iMpc~^), where the 
dependence is both radial and tangential. This is a prac- 
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Fig. 2. Slices in Rf{k) (solid, red) compared to R^{k) (dashed, black) for ^ = 5, 50 {top, bottom), for a wide and shallow 
galaxy survey (left, rp — 100/i~^Mpc, or ^ 0.05) and a wide and deep galaxy survey (right, ro — 1400/i^^Mpc or 
zm - 0.8). 



tical illustration of Equation 18 (see also Equation A. 5). 
This can be understood as tangential modes being atten- 
uated due to mode-cancelling along the line of sight - the 
attenuation becomes important when the selection function 
spans a wide redshift range. We refer to this effect as the 
radialisation of the power spectrum, where we have used 
the baryonic wiggles to probe the radialisation. We discuss 
the radialisation limit in more detail in §3.2 and in §4. 



One of the important results of Figure 1 is that the 
radial limit is still reached even when evolution (growth and 
bias) of the field are considered. This can only be considered 
numerically since there is no analytic solution to Equation 
31 due to the inclusion of the growth and bias terms. We 
test this again with strongly evolving values of the galaxy 
bias, e.g. b{z) — (Vl -|- z) , where n = —10,2, and find 
this still has little effect on the radialisation of Rf{k). The 
growth values for standard concordance cosmology do not 
seem to have an effect on reaching the radialisation limit. 



3.2. Radialisation Limit 

We arc interested in quantifying the effect of the radialisa- 
tion on the BAOs themselves. In Figure 2, we consider slices 
in £ and k through Rf (k) in order to investigate how rapidly 
the limit in Equation 18 is reached. The black (dashed) lines 
in Figure 2 correspond to R^{k) (which are independent of 
£ and of the survey selection function) and the red (solid) 
lines correspond to Rf{k) for different slices in £ space and 
for narrow (LHS) and wide (RHS) selection functions. As 
the survey selection function is widened, we find that: 

lim Rf{k) = R^{k), (34) 

i.e. Rf{k) tends towards the oscillations in R^{k), not only 
in phase, but also in amplitude. This is noticeable in Figure 
1 : the amplitude of the first peak for example is higher for 
the deep survey (> 9%) than for the shallow survey (< 8%). 
This is another check of Equation 7. 

Equation 19 gives us an analytic way to study the 
limit where we can expect the SFB power spectrum to be- 
have radially. In Figure 3, we plot the lines correspond- 
ing to the radial limit in {k,£) space for different surveys 
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become radial before the full SFB power spectrum does, es- 
pecially for large I. This surprising effect, further motivates 
the use of the SFB analysis for BAOs. 
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Fig. 3. Radialisation limit as defined by Equation 19 (di- 
agonal lines) for Gaussian surveys (Equation 15) with 
ro = 40, 70, 200, 300, 1400/i^iMpc, i.e. roughly correspond- 
ing to the following surveys: IRAS (^mean ^ 0.015), 2MRS 

V^mean 

0.03), SDSS (2„iean ~ 0.10), 2dF ( 0.15), 
and a stage IV survey (zmcan ~ 0.80). Information above 
the diagonal line is in the radial limit. We overplot the phys- 
ical scales corresponding to the first three BAO peaks for 
comparison. 

with Gaussian distributions (Equation 15) corresponding to 
ro = 40, 70, 200, 300, 1100 [/i-^Mpc] . We have chosen these 
values as they correspond roughly to the following sur- 
veys (for the chosen fiducial survey): IRAS (zmoan ^ 0.015, 
Fisher et al. 1995a; Strauss et al. 1992), 2MRS (2 MASS 
All-Sky Redshift Survey, z^ean ~ 0.03, Huchra et al. 2011), 
SDSS (Sloan Digital Sky Survey, Zmoan ^ 0.10, Adclman- 
McCarthy et al. 2008), 2dF (2 degree Field, z^oan 0.15 
Colless et al. 2003) and a Stage IV galaxy survey (zmcan ^ 
0.80, Albrccht et al. 2006). Note that, in principle, these 
surveys would have to be all-sky to obey the radialisation 
limits plotted in Figure 3. 

The area above each diagonal line corresponds the ra- 
dial limit, i.e. where the radialisation has taken place and 
Equation 18 holds. We overplot the radial scales corre- 
sponding to the first three peaks of the BAOs to show how 
these will be probed for different survey depths. 

By comparing the diagonal lines in Figure 3 with the 
BAO 'turnover' (i.e. the ^-dependent fc-scale at which the 
BAOs switch from being mostly tangential to mostly ra- 
dial) in Figure 1, we notice however that the BAOs seem to 



4. Conclusion 

In this paper, we have presented a new way to study Baryon 
Acoustic Oscillations (BAOs) using a Spherical Fourier- 
Bessel (SFB) decomposition of a wide-field deep galaxy 
survey. The BAO signal in SFB space can be studied by 
considering the ratio {k) of the SFB power spectrum 
with wiggles {C^{k)) to the smooth SFB power spectrum 
(C"°''(A:)), similarly to what is done in Fourier space. In 
this decomposition, BAOs can be observed simultaneously 
in both the radial [k) and tangential (£) modes. 

For a field which is statistically isotropic and homo- 
geneous (SIH), the SFB power spectrum is purely radial, 
i.e. independent of tangential modes, and reduces to the 
Cartesian Fourier power spectrum. In the other extreme 
limit where the field is in a thin shell, we showed that all 
the SFB information is contained in the tangential modes 
and is simply related to the 2D spherical harmonic power 
spectrum. 

Considering the practical case where the field is ob- 
served with a radial selection function (thus partially vi- 
olating the SIH condition), we find that the radial limit 
can still be reached for selection functions covering a large 
radial range (see Appendix A), which we refer to as the ra- 
dialisation of the power spectrum. The radialisation will 
be limited to a region in [t, k) space corresponding to 
k ^ ■y/^(£ + l)/r'o, where ro parameterises the radial se- 
lection function. We find that radialisation is a good ap- 
proximation for these modes even when evolution due to 
growth and bias are considered, both of which can be con- 
sidered as a further violation of the SIH condition. To study 
this limit, we have derived an analytic solution to the win- 
dow function of a non-evolving, unbiased galaxy field when 
the radial selection function is a Gaussian in r centered on 
the observer (Equation 16). 

We also find the BAOs considered in SFB space radialise 
as the survey depth is increased, meaning both the phase 
and amplitude of the BAOs tend towards the Fourier space 
ratio R^{k). This means that BAOs for a wide-field shallow 
survey have smaller amplitude and are spread across the 
(£, k) space, while BAOs for a wide-field deep survey have 
a larger amplitude and are confined to the radial modes 
(Figures 1 and 2). 

We study the radial limit analytically (Figure 3) and 
find that it can in principle be observed (for a limited 
range and for small physical scales) with all-sky surveys 
with current surveys depths or for future stage-IV surveys. 
This suggest that the first BAO peak becomes radial up to 
t ^ \{) for an all-sky survey with similar depth to SDSS 
or 2dF, and up to ^ ^ 70 for an all-sky stage IV survey. 
In practice though, the radialisation might be observable 
to even higher I since we observe that the BAOs radialise 
more rapidly than the full SFB power spectrum, especially 
at large Subsequent BAO peaks also become radial for 
large values of t and for shallower surveys, though these 
may be already be in the non-linear regime. For modes that 
have become radial, measurements at different .^'s are useful 
in practice, to reduce measurement errors due to cosmic 
variance and shot noise. 
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We note that we have ignored redshift-space distortions 
in our analysis though the prescription for these in SFB 
space is well known, (see for e.g., Heavens & Taylor 1995). 
These distortions may affect the radialisation as they will 
introduce mode-mixing, albeit with a distinct signature and 
are readily described in the SFB basis. Further ^-mode mix- 
ing will occur for incomplete sky coverage, which can be 
corrected for using the mask geometry. 
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Lanusse for help implementing this. We extended iCosmo" software 
for our calculations (Refregier et al. 2011). This research is in part 
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Appendix A: Radial Limit for a Gaussian Selection 
Function 

As shown in §2.2.1, when the galaxy selection is Gaussian, 
i.e. (j){r) = e"^''/''"-' , it is the galaxy selection function 
(Equation 13) reduces to 



Weik,k') = 



■ exp 



]kk' 



(A.l) 

where /^(x) is the modified Bessel function of the first kind. 
For numerical reasons, it is often useful to evaluate the 
exponentially scaled modified Bessel function of the first 
kind I^ix) = exp(— a;)/iy(a;) instead of I^{x), in which case 
it is useful to re-write the above Equation as: 



Weik,k') = 




exp 



k~k' 



\kk' 



(A.2) 

Figure A.l shows the window function for two values of rg, 
£ = 3 as a function of fc, k' . As rg becomes large, it can 
be seen that the window function tends towards a delta 
function -p5{k — k'). 

To study this limit more precisely, we use the asymp- 
totic form for the modified Bessel function I„{x] ~ 



X ^ — 1/4|, which gives 



WfXk.k') 



\/7r r-Q 



k-k' 



(A.3) 



in the limit where r^kk' ^ 2i{l +1), or approximately. 



rok > ^J2^{l+l). 



(A.4) 



Using the definition of a Dirac delta function, namely 
J dx h{x)5{x — xq) = h(xo) for an arbitrary function h{x), 
the window function becomes, in the limit j'g — !■ oo. 



We{k,k')c,^6{k-k'), 



(A.5) 



in agreement with Equation 17 for the radial case discussed 
in §2.2.2. Thus, if the condition of Equation A.4 holds, the 
power spectrum Ci^'^{k, k') becomes radial (i.e independent 
of t) and equal to Ct{k)5{k - k!) = P{k)S{k - k') as in 
Equation 18. 
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A. Rassat ^'^ and A. Refregier ^: 3D Spherical Analysis of BAOs 




Fig. A.l. Window function Wg{k, k') for the Gaussian selection function with for £ = i and tq = lOO/i ^Mpc (left) and 
To — 1400ft-^^Mpc (right). As the selection function parameter tq becomes larger (i.e., increasing the redshift coverage 
from Znicd ^ 0.05 to ^mcd ^ 0.8), the window function tends towards ■^S{k' ~ k) for scales where k ^ y^'2,£{£ + 1) /ro- 
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